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The chiral properties of lattice fermions can be improved by altering either their fermion-gauge 
coupling or the pure gauge part of the action (or both). Using both perturbation theory and 
nonperturbative simulation, we compare a simple alteration of the gauge action (which encompasses 
the Wilson, Symanzik, Iwasaki, and DBW2 actions), and HYP-blocked links in the fermion action. 
Perturbative tests include calculations of the potential, flavor-changing quark scattering amplitudes, 
and matching factors for currents. Non-perturbative tests include the potential, measurements of 
flavor symmetry breaking for staggered fermions, the behavior of topological objects, and properties 
of overlap actions. Our results display the bad properties of these actions as well as their good ones. 


I. INTRODUCTION 

The chiral symmetry breaking of lattice fermions is the most serious lattice artifact we face in numerical simulations. 
Most dynamical simulations use staggered or Wilson-like fermions, where chiral symmetry breaking is explicit and 
difficult to control. In the case of staggered fermions lattice artifacts induce flavor symmetry breaking, resulting in 
non-Goldstone pions that can be significantly heavier than the Goldstone particle. In the case of Wilson-like fermions 
the effect of chiral symmetry breaking is most observable in the occurence of exceptional configurations. Domain wall 
fermions, which in theory are chiral, have a residual chiral symmetry breaking due to the finiteness of the fifth lattice 
direction. Overlap fermions are chiral, no matter what non-chiral action they are constructed from. However, the 
computational requirements needed to evaluate the overlap operator can be significantly lowered if it is based on a 
near-chiral action. 

A systematic improvement program modifies the pure gauge part of the action, the fermionic action, and the gauge 
connections of the fermions, to construct a chirally improved action with reduced lattice artifacts. Less ambitious 
programs might change only part of the action at a time, according to the improvement needs of the program. 

Several possibilities have been explored to improve the chiral symmetry of fermionic actions. They are all based 
on the observation that most of the chiral violation come from the short scale, plaquette level vacuum fluctuations of 
the gauge field. One option is to smear the gauge connection of the fermions, reducing the coupling of fluctuations 
to the fermions. This option has been used with staggered fermions i I- The two smearing transformations that 
are used most extensively are the Fat7/Asqtad jH ^ and HYP smearings |^. Both of these transformations are 
O(a^) perturbatively improved, though only the Asqtad smearing was designed perturbatively. The HYP smearing 
is optimized non-perturbatively. 

The other option that has been studied lately is to modify the pure gauge action such a way that the creation 
of small scale vacuum fluctuations is dynamically suppressed. The Iwasaki action [Q has been used extensively by 
CP-PACS both in dynamical Wilson and domain wall fermion simulations H, while the DBW2 action @,0 was 
chosen by the Columbia-BNL group in their domain wall fermion simulations ||^. The Iwasaki and especially the 
DBW2 actions have larger perturbative lattice corrections than the Wilson plaquette action though in numerical tests 
that did not seem to increase the lattice artifacts. 

The advantage for numerical simulations of modifying the gauge action instead of smearing the gauge connection 
is obvious. Almost any gauge action can be simulated faster than a complicated fermionic action, though with the 
recently developed partial-global stochastic Metropolis (PGSM) update even projected smeared link fermions can be 
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simulated effectively I, SOI- The important question is if the improvement offered by the modified gauge actions is 
sufficient, and even more, if, in addition to the improved chiral symmetry any unwanted lattice artifacts are introduced 
by these actions. 

Our goal in this paper is to compare different physical properties, both perturbative and non-perturbative, of the 
smeared and modified gauge actions, and possibly to predict which choice is going to give the most efficient approach 
of the continuum limit. 

In Sec. 2 we introduce the actions we study. A perturbative analysis of the properties of these actions is carried 
out in Sec. 3. In Sec. 4 we show simulation results for the heavy quark potential. The topological properties of these 
actions are discussed in Sec. 5. Flavor symmetry violations for staggered actions are shown in Sec. 6. Sec. 7 discusses 
some properties of overlap actions built using HYP links or in the background of gauge fields with the gauge actions 
we study. Our conclusions are found in Sec. 8. 


II. THE ACTIONS 


In the following we consider a family of actions consisting of the 1x1 plaquette and the 1x2 planar loop 

Sg{U) = ^[co Y. + E ( 1 ) 

with the normalization condition co + 8ci = 1. The coefficient ci can vary, giving the different actions 

Cl = 0 Wilson (2) 

— 1/12 tree level Symanzik (3) 

—0.331 Iwasaki (4) 

-1.4088 DBW2. (5) 

The tree level Symanzik action is 0[a^) improved 0, while both the Iwasaki and DBW2 actions have opposite 0{a?) 
corrections than the Wilson plaquette action. Some of the scaling properties of the Iwasaki action have been studied 
in 1^ while the DBW2 action was investigated in 0|ii- In our numerical symulations we use the 1-loop tadpole 
improved Symanzik action .0 
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where is a six-link parallelogram with links running around the opposing edges of the cube. The /3 coefficients 
for S'symii are tadpole improved according to the plaquette expectation value. 

In addition to the different gauge actions we will consider two smearing transformations for the fermion-gauge field 
interaction. The Fat7 smearing replaces the usual one-link coupling with a linear combination of gauge loops up to 
length seven. It has perturbatively determined coefficients which remove the flavor changing gluons at the edges of 
the Brillouin zone for staggered fermions. When an additional 5-link term is added to the Fat7 smearing the resulting 
action is O(a^) improved in the fermion-gluon connection ||^. With tadpole boosted coefficients this leads to the 
Asqtad smearing transformation ||^. (One should note that the Asqtad staggered action also has a third nearest 
neighbor Naik term, which we do not include here.) The Asqtad smeared links are not unitary - they are simply the 
linear combinations of the extended gauge paths. 

Our second smearing transformation is HYP smearing |^. HYP smeared links are constructed from three levels 
of modified, SU(3) projected APE blocking steps in a way that makes the transformation local and smooth. While 
the HYP smearing is non-perturbatively optimized, its coefficients can just as well be tuned to give perturbative 
improvement. 





3 


III. PERTURBATIVE CONSIDERATIONS 
A. Preliminaries 


The gauge actions we study have propagators which obey the equation 


^ kf^6p^Qfj,pkp')^Gfj^jy — Spj^ 


p 


( 7 ) 


where 


= (1 - - ci{kl + kl)) 


( 8 ) 


and kp = 2sin(fc^/2). ^ is the gauge-fixing term: ^ 


1 is Feynman gauge. The Wilson gauge action propagator is 


G,. = [<5^. + (?-l)^]^. (9) 

For other actions we merely numerically invert Eq. to construct the propagator. 

We will be concerned only with unitary fat links, gauge connections which are themselves elements of the gauge 
group, even though they may be built of sums of products of the original thin links of the simulation. For smooth 
fields the fat links have an expansion Vp{x) = 1-1- iaBp{x) -I- ... and the original thin links have an expansion 
Up(x) = 1-1- iaAp(x) -I-.... For computations of 2- and 4-quark operator renormalization/matching constants at one 
loop, only the linear part of the relation between fat and thin links is needed [^, and it can be parameterized as 

^ hp^{y)A^{x + y) . (10) 

Quadratic terms in (|l^ , which would only be relevant for tadpole graphs, appear as commutators and therefore do 
not contribute, since tadpoles are symmetric in the two gluons. In momentum space, the convolution of Eq. ( 0 ) 
becomes a form factor 

BM) = hp^{q)A^{q) . ( 11 ) 


If all gluon lines start and end on fermion lines, then, effectively, the gluon propagator changes into 

hpxG\a h(ji /. ( 12 ) 


Obviously, this means that in perturbation theory, as far as the fermions are concerned, fattening the links in the 
fermion action is equivalent to altering the gauge action. 

The family of smearings including HYP links and the order-a^ improved link 1^ which, when augmented by 
tadpole improvement, gives the “Asqtad” link used by the MILC collaboration ||^, have 

hpu = 5y^pDp{k) -I- (1 - 5v,p)Gy^p{k). (13) 


The diagonal and off-diagonal couplings can be decomposed, respectively, as 
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where all indices {p,,v,p,<j) are different. 

The coefficients di _4 distinguish the different choices of links: 


1. Fat-7 links: 


2. 0{A) improved links: 


di — 1 , (^2 — 1 , d^ — 1 , d/^ — 0 . 


di — 0, c?2 — 1, da — 1, ^4 — 1- 


(15) 

(16) 


(17) 


(18) 


3. HYP smeared links: 

di = ( 2 / 3 )q!i (1 -|- 0:2(1 + 03)), d2 = (4/3)0102(1 -f 203), da = 801O2O3, d4 = 0 . ( 19 ) 

There are two interesting choices for the at. The first was determined in Ref. using a non-perturbative 
optimization procedure: Oi = 0.75, 02 = 0.6 03 = 0.3. The second is chosen so to remove 0{A) flavor- 
symmetry breaking couplings at tree level. This gives the same di as for Fat-7 links. 


B. Static Potential 

With the gluon propagator we can immediately compute the static Coulomb potential, 

Vc{r) = f exp{ik ■ r)Goo{qfj. = {0,k)) (20) 

J ak 


where = Ylj d(o.kj)/(2Tr) will be the symbol for integration over the (rescaled) momentum hypercube. With 

our conventions, the continuum potential is V(r) = l/(47rr), and so plotting the rescaled lattice potential 47rrP(r) 
immediately exposes the lattice artifacts of a particular action. We show results for this quantity in Fig. |]for the four 
candidate actions of our study. The tree level Symanzik action has no 0(A) nor 0(A} discretization errors and also 
has the smallest scaling violations. The other actions have 0(A) scaling violations, which should be of opposite sign 
for the Wilson action versus the Iwasaki and DWB2 actions. This is seen in the figure. However, the most noticable 
feature of the potential is the systematic distortion of iirrV (r) below unity at small lattice distance as the coefficient 
Cl becomes more negative. The results of this calculation strongly disfavor use of a large negative value of Ci for any 
lattice simulations with physics related to the short distance part of the potential. 

The potential itself is, of course, unaffected by any fattening of the fermion’s gauge connection, but we can also 
define a “smeared potential,” in which the gluon propagator Gqo is replaced by Qoo- This quantity has the physical 
interpretation of the potential seen by a heavy lattice quark whose gauge connection is a fat link. Results for two 
gauge actions (Wilson and tree level Symanzik) and two smearings, HYP and Asqtad, are shown in Fig. 

Both of these smearings distort the lattice potential at small r/a. The immediate conclusion that one draws from 
these pictures is that one should not do simulations involving heavy quarks with smeared links - the loss of Coulomb 





5 






Figure 1: The scaled static potential 47rrV^(r) for several lattice actions: (a) Wilson (ci = 0), (b) tree level Symanzik (ci = 
-1/12), (c) Iwasaki (ci = -0.331)(d) DWB2 (ci = -1.4088). 


behavior at short distance will allow the wave function of a heavy quark-antiquark bound state to spread out, the 
value of the wave function at the origin will be small, and hyperfine splittings will collapse. This effect is readily seen 
in simulations. 

A particular example of the danger of fattening heavy quarks is seen in the recent MILC study of heavy quark-light 
quark decay constants, fo and /b ||1^ . One of the data sets collected by these authors used fermions fattened with 
a large amount of APE-smearing |20|| (c = 0.45, Al = 10 in the conventions of Ref. |j^). This amount of smearing 
produces a noticable suppression of the static potential out to r/a ~ 3 — 4. The authors observed a twenty per cent 
reduction in Jb- This effect presumably would go away at smaller lattice spacing, but comparing Fig. ^ one would 
need to halve the lattice spacing to reduce the lattice artifact to the level of HYP blocking. 

However, the motivation for using smeared links in a fermion action is, by and large, to improve the chiral properties 
of the action. This is physics applicable to light quarks, not heavy ones. Many simulations of light quarks with various 
degrees of fattening show no ill effects on spectroscopy or on matrix elements - generally, improvement of scaling is 
observed ||^, . Fattening is something which can be done selectively, in a mass-dependent way. Altering the gauge 

action will affect quarks regardless of their mass. 

Finally, we call the reader’s attention to the rather large lattice spacing artifacts of the Asqtad-smeared potential 
as compared to the HYP- smeared potential (or any of the usual unsmeared potentials). 
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Figure 2: The scaled static “smeared potential” 47rr’F(r) for two lattice actions and two smearing functions: (a) Wilson (ci = 0) 
gauge action with HYP blocking, (b) tree level Symanzik (ci = —1/12) gauge action with HYP blocking, (c) Wilson (ci = 0) 
gauge action with Asqtad blocking, (d) tree level Symanzik (ci = —1/12) gauge action with Asqtad blocking, 


C. One-loop Perturbative Matching Factors 

A matrix element of an observable computed using one regularization (MS) is related in perturbation theory to 
a linear combination of observables computed using another (lattice) regularization through a matrix of matching 
coefficients Znm 


(/|Or‘(M)l*)M5 = E a)„™(/|0'““(a)™|t) (21) 

m 

where each coefficient is a difference between (the finite part of a) continuum-regulated and a lattice-regulated ex¬ 
pression 

= 1 + ( 22 ) 

is the squared coupling). For currents, it is customary to divide out the quadratic Casimir Cp and to present 
results for z = ~ ^iatt)/CF- For four-quark operators (for weak-interaction matrix elements, for example), 

it is customary just to write b — A|^ — Aiatt- Generally, in the context of perturbation theory, one attempts to 
design actions, operators, and methods of performing perturbative calculations so that the z or b coefficients are 
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Wilson 

Tree level Sym 

Iwasaki 

DWB2 

Zv 

-15.33 

-11.91 

-7.44 

-3.03 

Za 

-13.79 

-10.72 

-6.71 

-2.75 

Zs 

-19.31 

-15/08 

-8.90 

-1.04 

Zp 

-22.38 

-17.47 

-10.36 

-1.59 

z+ 

-36.63 

-28.86 

-19.15 

-1.08 

Z- 

-43.20 

-32.78 

-18.27 

-1.61 


Table I: a and b coefficients for one loop matching factors for some two- and four-fermion operators, for thin-link clover fermions 
with Csw = 1, and gauge action as labelled. Errors are ±1 in the last digit shown. 



Wilson 

Tree level Sym 

Iwasaki 

DWB2 

Zy, Za 

-14.80 

-11.24 

-6.76 

-2.63 

Zp, Zs 

-39.24 

-30.35 

-17.99 

-4.38 

z+ 

-25.18 

-19.25 

-12.54 

-7.86 

Z- 

-68.06 

-51.45 

-29.0 

-5.35 


Table II: 2 and b coefficients for one loop matching factors for some two- and four-fermion operators, for naive fermions with 
Csw = 1, and gauge action as labelled. Errors are ±1 in the last digit shown. 


minimized. The reader might recall that at typical lattice spacings studied today, typical choices for g^/47r range 
from 0.1-0.2, so a z or 6 of about 20 implies a twenty to forty per cent effect from a one loop calculation-perhaps a 
bit large for comfort. Typically, with fat links one can reduce these numbers by an order of magnitude, as seen by 
perturbative calculations |^, or nonperturbative simulation |Q. 

We have computed one loop perturbative matching factors for local currents 'ip{x)T'ip{x)-ior the vector (V), axial 
vector (A), scalar (S) and pseudoscalar (P) currents, as well as the matching factors for the four-fermion operators 
0± = Oi ± O 2 , built from 


o = ) 0 (gt)r2gf). 


(23) 


where Pi = r2 = 7^(1 — 75)- 0 ± = Oi ± O2', if color labels a = d, P = j, O = Oi; if color labels a = P, 'j = 6 , 
O = 02 - (These operators have no penguin contributions.) 

An important indicator of chiral improvement ||^ for nonchiral actions is the difference zy — za = {zp — zs)l‘2.-. 
the mixing of four-fermion operators into the opposite chirality sector is controlled by this quantity. 

Our results are shown in Tables li For operators with anomalous dimensions (all but the vector and axial vector 
currents) our results are for the case (lattice spacing ax regularization point /i) = 1. We see that as ci becomes 
increasingly negative, the z and b coefficients generally shrink. 

We next repeat these calculations, but now with fermions with HYP-smeared gauge connections. Tables HI, IV 


show a dramatic reduction in the z and b coefficients, even when the Wilson gauge action is used. With the reduction 
of these numbers comes also a reduction in the difference zy — z^ for clover fermions. 

From the point of perturbative theory for matching factors, the conclusion is clear: it is much more efficient to 
fatten the fermion gauge connections than to increase ci in the gauge action. 
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Wilson 

Tree level Sym 

Iwasaki 

DWB2 

Zv 

-1.38 

-1.18 

-0.89 

-0.50 

Za 

-1.30 

-1.11 

-0.84 

-0.48 

Zp 

0.04 

0.54 

1.55 

3.77 

Zs 

-0.12 

0.41 

1.45 

3.72 

z+ 

-6.43 

-6.14 

-5.89 

-6.12 

Z- 

2.16 

3.11 

4.84 

8.34 


Table III: z and b coefficients for one loop matching factors for some two- and four-fermion operators, for clover fermions with 
Csw ~ 1, HYP-smeared links, and gauge action as labelled. Errors are ±1 in the last digit shown. 



Wilson 

Tree level Sym 

Iwasaki 

DWB2 

Zv, Za 

-0.95 

-0.79 

-0.59 

-0.33 

Zp, Zs 

-0.62 

-0.01 

1.20 

3.61 

z+ 

-4.75 

-4.65 

-4.75 

-5.53 

Z- 

1.92 

2.95 

4.81 

8.37 


Table IV: z and b coefficients for one loop matching factors for some two- and four-fermion operators, for naive fermions with 
Csw = 1, HYP-smeared links, and gauge action as labelled. Errors are ±1 in the last digit shown. 


D. On-shell Scattering Amplitudes 

Let us recall that the motivation for introducing smeared links into staggered fermions was to suppress the coupling 
between the region of the fermion Brillouin zone near = (0, 0,0,0) and regions of the Brillouin zone corresponding 
to doublers: one or more kfj. ~ tt. Smearing amounts to a form factor which suppresses the emission or absorption 
of gluons whose exchange could scatter a quark from a ~ 0 into a doubler state. The absence of this kind of 
scattering means an improvement in flavor symmetry for staggered fermions since a process which transforms a quark 
of one flavor (living in one part of the Brillouin zone) into another flavor is reduced. 

To quantify this scenario, let’s imagine (in the continuum, first) the scattering of two on-shell quarks of momentum 
±pi into two quarks of momentum ±p 2 - The T—matrix for the scattering is 

T{k) = [u{p2)-f^,u{pi)]G^i,{k)[ui-p2h„ui-pi)] (24) 

where of course k = pi — p 2 - Imposing the on-shell constraint 7 • p\u{pi) = 0, we see that the gauge-dependent term 
(proportional to (^— l)A:^fc^/fc^) vanishes, leaving the scattering amplitude proportional to the Feynman gauge gluon 
propagator (and spinor factors). The reader can quickly confirm that the same result obtains for naive fermions and 
the Wilson gluon propagator. 

Now recall that the effect of smearing is to replace the gluon propagator by the “smeared gluon propagator,” Eq. 
[ 1 ^ . As a way of comparing the level of flavor symmetry violation, we just look at the Feynman gauge propagator 
(appropriate for flavor symmetry restoration with thin link fermions) or the smeared propagator. We plot k^ 

versus for a set of momenta = nn^/S, = 1.. . 8 . Fig. || shows our results for thin link actions. As ci 
becomes more negative, the “scattering amplitude” at large k decreases. We would anticipate that flavor symmetry 
would be improved by negative ci. We will confirm this expectation in Sec. 6 . 

However, much greater suppression at large k is achieved by converting to a smeared link. In Fig. ^ we show several 
examples of G)j.)j., which parameterizes the scattering of smeared link fermions. HYP links and the Wilson gauge 

action produce more suppression than thin links and the DWB2 action. Clearly, a combination of greater ci and a 
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Figure 3: fe^G^^(fc) for several lattice actions: (a) Wilson (ci = 0), (b) tree level Symanzik (ci 
(ci = -0.331)(d) DWB2 (ci = -1.4088). 


— 1/12), (c) Iwasaki 


smeared link would produce a larger effect. 

There is an obvious qualitative connection between Figs. ^ and Q and the results of one-loop perturbation theory 
for matching coefficients: As the magnitude of the gluon propagator shrinks at large fc, so does its contribution to the 
integrals of the perturbative calculation. 


IV. THE NON-PERTURBATIVE STATIC POTENTIAL 

We start our discussion of the non-perturbative properties of the different actions with the static quark potential. 
The perturbative results of Sect. 3.2 suggest that the thin link Symanzik action has the smallest lattice distortion. 
The Wilson action has a positive lattice correction at small distance and observable rotational symmetry breaking 
even at rja = 2 — 3. The perturbative DBW2 potential has a larger and negative correction at rja = 1 and observable 
rotational symmetry breaking even at rja = 4 — 5. The perturbative HYP smeared Wilson gauge action potential 
has similar distortion at rja = 1 as the DBW2 but much smaller rotational symmetry violation while the Asqtad 
smeared Wilson gauge potential has smaller distortion at short distances but larger rotational symmetry violation. 

We studied the non-perturbative static potential on Wilson, 1-loop Symanzik and DBW2 gauge backgrounds with 
and without HYP smearing. The simulations were carried out on 8^ x 24 lattices with Sommer scale r^ja ^ 3.0 for 
all three actions, at /3 = 5.7 for Wilson, f3 = 7.775 for Symanzik and (3 = 0.82426 for the DBW2 action. At large 
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Figure 4: k^Q^p,{k) for two lattice actions and two smearing functions: (a) Wilson (ci = 0) gauge action with HYP blocking, 
(b) tree level Symanzik (ci = —1/12) gauge action with HYP blocking, (c) Wilson (ci = 0) gauge action with Asqtad blocking, 
(d) tree level Symanzik (ci = —1/12) gauge action with Asqtad blocking. 


distances the smeared and thin link potentials differ only by an irrelevant additive constant. In Fig. ^ we plot both the 
smeared and thin link potentials, the former one shifted by a constant to match the thin link potential at r/a = -s/T, 
a somewhat arbitrarily chosen matching point. We fit each potential following the method proposed in and used 
with the HYP potential in ||^, taking a four parameter functional form 

V\s,tt{r) = V/ont(c) + eiydr) - -^) 

47rr 

where V/ont is the continuum potential 

e 

Vcont =-h Vo + err, 

r 

and Vc{r) is the lattice Coulomb potential of Eq. ^ The term e(Vc(r) — l/(47rr)) is an attempt to model and remove 
the lattice artifacts of the potential. It is difficult to judge how much of the lattice artifacts can be described by this 
term, only the quality of the fit can justify its use. 

Fig. ||a shows the Wilson gauge action potential measured with thin links (diamonds) and with HYP links (oc¬ 
tagons). The dotted line of the ffgure correspond to the fit of the continuum potential V/ont after the removal of 
the lattice artifacts while the dashed line is the same continuum potential obtained with the HYP smeared links 
and their corresponding perturbative corrections. The fact that it is impossible to resolve the two different lines in 
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Figure 5: The static potential measured with a) Wilson gauge action, b) 1-loop Symanzik gauge action and c) DBW2 gauge 
action. In all cases both the thin link (diamonds) and HYP smeared (octagons) potentials are plotted, shifted to agree at 
r ja — V7. The dotted and dashed lines are the fitted continuum potentials as described in the text. 


the figure indicates that the lattice artifacts are consistently removed. Also, the sign and relative magnitude of the 
lattice corrections are what we expected from the perturbative Coulomb potential. Fig. lb is the same for the 1-loop 
Symanzik action. The smeared and thin link results are consistent, indistinguishable. The thin link potential has 
very small lattice artifacts but after removing the lattice correction both thin and HYP potentials predict the same 
continuum values. In Fig. ||c we plot the corresponding potential data for the DBW2 action. The agreement between 
the thin and HYP smeared potentials is good though not as perfect as in the previous two cases, mainly because of the 
stronger rotational symmetry violation of the thin link potential.The results agree with the perturbative predictions: 
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Figure 6: The action of smooth instantons, normalized by the continuum value, as the function of the instanton radius, 
calculated with different actions, a) solid line: Wilson action, dotted line: tree level Symanzik action; dashed-dotted line: 
Iwasaki action; dashed line: DBW2 action, b) solid line: Wilson action (thin link); dashed line: HYP smeared Wilson ation; 
dotted line: HYP smeared tree level Symanzik action. 


even the thin link DBW2 potential has a large distortion at r/a = 1, about the same as the HYP link potential with 
Wilson or 1-loop Symanzik action. If the value of the potential (or any other quantity) is important at rja = 1, the 
DBW2 action is not a good choice to use. 

It is difficult, if not impossible, to prove how the lattice Coulomb term could describe the lattice artifacts of the 
non-perturbative potential. Only the fact that the different potential measurements with thin and HYP links give 
consistent continuum results for all three gauge actions justifies its use. 

Recent calculations of the static potential with the DBW2 action did not show such a large distortion at small 
distances. However in Ref. only the on-axis potential was measured at distances r/a > 2. From that data it 
would have been hard to see the distortion at small distances. 


V. TOPOLOGY AND THE GAUGE ACTION 
A. Smooth instantons 

It is generally believed that flavor symmetry violation for staggered fermions, the residual chiral symmetry violation 
of domain wall fermions, and the computationally most demanding small eigenmodes of the overlap fermions are due 
to small scale, plaquette level vacuum fluctuations. Smearing attempts to remove these objects by averaging the 
gauge links while the Iwasaki and DBW2 actions do the same by making it energetically unfavorable to create them 
in the first place. That mechanism can be seen clearly from the action of smooth instantons. We calculated the 
instanton action on a set of smooth instantons with varying radii. These instantons were created in singular gauge 
on 32"* lattices and blocked twice in order to approximate the smooth, continuum solution. In Fig. ||a we show the 
result, normalized by the continuum instanton action, for Wilson, tree level Symanzik, Iwasaki and DBW2 actions. A 
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Figure 7: The Monte Carlo time history of the topological charge with various gauge actions. The actions are (from top to 
bottom) Wilson (/3 = 6.0), Symanzik (/3 = 8.4), Iwasaki {(3 = 2.6) and DBW2 (/3 = 1.04). 


perfect lattice action should have a profile that is less than one for small radii and one for radii larger than a critical 
value where the vacuum fluctuation becomes an instanton. At the critical radius the action function is non-analytic, 
it develops a kink. None of our actions can reproduce this behavior, but not surprisingly the closest to it is the tree 
level Symanzik action. The Wilson action profile approaches the continuum value more slowly, from below. The 
Iwasaki action overshoots the continuum value by about 50% at rja^ 0.6, suppressing fluctuations of this size. The 
correponding DBW2 curve is nothing less than shocking. The curve rises to almost four at rja ^ 0.6 and even at 
distance rja = 1.5 it is above two. The DBW2 action strongly suppresses instantons and dislocations with radius 
0.3 < rja <2 — 3. The very small fluctuations are still present, but small radius instantons are disfavored. On lattices 
where these small instantons are important physically one would expect fairly large lattice artifacts from the DBW2 
action. 

Smearing attempts to remove dislocations seen by the fermions by averaging the gauge links. Fig. compares the 
HYP smeared instanton profiles of the Wilson and Symanzik actions. For reference we include the thin link Wilson 
action profile again. The most important conclusion from Fig. ^ is that smearing removes dislocations with radius 
rja < 0.5, the action profile rises only at r/a ~ 0.7. Even though the gauge configuration can have plenty of small 
instantons and dislocations, most of these are not seen by the fermions. 


B. The autocorrelation of the topological charge 


We have seen that as the coefficient ci in the action becomes more negative, the action favors small instantons less 
and less. In a Monte Carlo simulation with local (one-link) updates, change of topology always occurs through the 
(dis)appearance of small topological objects. It is thus not very surprising that the suppression of small topological 
objects implies that topology changes less often. In Fig. we show the Monte Carlo time history of the topological 
charge with different gauge actions. The charge was measured using the RG improved charge operator ^ after 
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P P 

Figure 8: Instanton size distribution of the different actions after two HYP smearing steps, a) Wilson (dotted lines) and Iwasaki 
(dashed lines) actions; b) Wilson (dotted lines) and DBW2 (solid lines) actions. 


8 levels of APE smearing steps. The units on the horizontal axis correspond to ten full sweeps of a combination of 
one overrelaxation and one Metropolis step over the whole lattice. The lattice size in these simulations was 12^ in 
all cases and the [3 values were matched to correspond to the same lattice spacing, a = 0.095 fm, set by the Sommer 
scale. Indeed, the difference among the gauge actions is striking. Notice that for the DBW2 action the MC time 
scale is an order of magnitude different. The integrated autocorrelation time of the topological charge was estimated 
to be 100 and 700 sweeps for the Wilson and the Iwasaki action respectively. In the case of the DBW2 action the 
autocorrelation time is so enormous that the available data was not enough even to estimate it. It is also interesting 
to observe that the plaquette autocorrelation time is 10, 7 and 5.5 for the Wilson, Iwasaki and the DBW2 action. The 
very large autocorrelation time of the topological charge of the DBW2 action was also noted by the RBC collaboration 


It is also noticeable that as ci becomes more negative, and the change between topological sectors occurs less 
frequently, the charge also becomes closer on the average to integer values. This also indicates the strong suppression 
of small instantons and in general the suppression of gauge configurations close to the boundary between two charge 
sectors. 


C. The instanton size distribution 

The instanton action profiles of Fig. indicate a slight suppression of small instantons for the Iwasaki action, and 
strong suppression of small and even larger instantons for the DBW2 action. On configurations with a ~ 0.095 fm the 
average instanton radius is about r/a = 3 but smaller instantons should be also present. In order to see if the different 
gauge actions have different instanton size distributions we have measured instanton sizes on a set of a ~ 0.095 fm 
configurations. We measured the topological charge density after 2-4 levels of HYP smearing and compared it 
to smooth instanton profiles. This is the same method we used in Ref. 0 0 except we do not extrapolate the 
instanton size to zero smearing level. In Fig. ^ we compare the instanton size distribution of the different actions after 
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two HYP smearing steps. We have used the same 12^, a = 0.095 fm configurations we analyzed in the previous section 
and normalized the distribution by the number of configurations. As Fig. ^ illustrates, there is not much difference 
between the Wilson and Iwasaki actions. The size distribution peaks around p/a = 3.5 and both actions predict the 
same topological density. (We are not concerned about the physical significance of the topological density here. Since 
we have used the same lattice spacing and analysis method with both configuration sets, comparing the two densities 
gives information about the two actions.) In contrast, the Wilson and DBW2 actions differ significantly, as Fig. ||b 
illustrates. The smaller instantons are suppressed by the DBW2 action, and the topological density is about 30% less 
on the DBW2 than on the Wilson configurations. If we consider the physical picture of quark propagation, where the 
quarks hop from instanton to anti-instanton in the vacuum ||3^, lack of instantons could point to observable scaling 
violations in the light hadron spectrum. 


VI. FLAVOR SYMMETRY VIOLATION IN STAGGERED ACTIONS 

Smeared actions are used with staggered fermions because they considerably reduce flavor symmetry violations. 
Both the Asqtad and HYP smearing are O(a^) perturbative improved though the coefficients of the HYP smearing are 
non-perturbatively optimized. Relative to the thin link staggered action, Asqtad fermions improve flavor symmetry 
by a factor of five, HYP fermions by about a factor of ten. Based on our perturbative and instanton analyzes, flavor 
symmetry could also be improved by modifying the gauge action. At first this approach might look attractive: it is 
much easier to simulate thin link fermions with a complicated gauge action than smeared link fermions. However taking 
the easy way might have serious consequences later on. The potential data indicates that the DBW2 action distorts 
short distance behavior, the time evolution of the topological charge points to unacceptably long autocorrelation 
times, the low topological density could imply large lattice artifacts. Nevertheless, in this section we consider the 
possibility of using different gauge actions with and without smearing in staggered fermion simulations and investigate 
the level of flavor symmetry violation these actions show. 

We have studied the quenched staggered spectrum on our 8^ x 24, ro/a ^ 3.0, a ^ 0.17 fm configurations. Before 
presenting our results for the spectrum, first we look at the distribution of the plaquette on these configurations. In 
1^ it was argued and illustrated that the end tail of the plaquette distribution is correlated with flavor symmetry 
breaking of staggered fermions. The argument is quite simple: flavor symmetry violation is caused by the strongly 
fluctuating gauge fields at the hypercubic level. These gauge links create plaquettes with very small value, therefore 
the number of plaquettes with very small value indicate the level of flavor symmetry breaking the fermions observe. 

In Fig. H we show the tail of the plaquette distribution, normalized by the number of configurations for the Wilson, 
1-loop Symanzik and DBW2 gauge actions. Fig. ^a compares the plaquettes constructed from thin links. Not 
surprisingly, the DBW2 action is a factor of eight better than the Wilson action. The 1-loop Symanzik action is also 
better than the Wilson action by about a factor of two. In Fig. ^ we plot the plaquette distribution of the Wilson 
and 1-loop Symanzik actions after one level of HYP blocking. The latter is again a factor of two better than the 
former, but both are an order of magnitude better than the thin link DBW2 action. (Observe the scale difference of 
the two figures.) Does the implication from the tail of the plaquette agree with the actual flavor symmetry violation 
of the different actions? We use the parameter 


A 


TT 


- mg 
mo 


(25) 


which measures the relative difference between the Goldstone pion mass and the non-Goldstone pions, to compare the 
different actions. This quantity diverges at zero quark masses and depends strongly on the lattice spacing, but since 
we have done all the simulations at approximately identical lattice spacings and volumes, as a function of the 
Goldstone particle, is a good indicator of the flavor symmetry violation of the different actions. The results, shown 
in Fig. 1^, confirms these expectations. The thin link 1-loop Symanzik action improves flavor symmetry relative to 
the Wilson action by about 30%. The thin DBW2 action is even better, it has flavor symmetry violation at the level 
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Figure 9: The tail of the plaquette distribution for different gauge actions a) with thin links and b) after HYP smearing. Dotted 
line: Wilson, dashed line: 1-loop Symanzik, solid line: DBW2 action. Observe the scale difference of the two figures. 


of the Asqtad smeared Wilson action. HYP smearing, even on Wilson configurations, is a factor of two better. One 
can reduce flavor symmetry violation even further by using the 1-loop Symanzik action. It is not easy to see from the 
figure, but 1-loop Symanzik is about 30% better than Wilson, even after HYP smearing. Since the Symanzik action 
has better scaling scaling properties than the Wilson action and does not suffer from topological autocorrelation slow 
down like the DBW2 and Iwasaki actions do, a HYP smeared Symanzik action appears to be the best choice for 
staggered simulations. 


VII. SMEARING, IMPROVED ACTIONS AND THE OVERLAP OPERATOR 

A. Fermionic charge and overlap 

Since a large negative ci suppresses the creation of small instantons or other objects where the gauge field is on 
the boundary between different topological sectors, gauge actions with a more negative ci result in smaller residual 
masses in domain wall simulations. A similar mechanism is at work in the case of the overlap. The overlap operator 
is defined in terms of a simple Dirac operator (e.g. Wilson) Dq by the formula 

Dov = I - A[A'^ AY^ , A=l + s-Do, (26) 

where s is a real parameter. If the gauge configuration is close to the boundary between different topological sectors, 
A^A has to have a small eigenvalue. Therefore, the suppression of these boundary gauge configurations can also thin 
out the small eigenvalues of A^A. This is important for practical applications since the cost of the overlap is governed 
by the condition number of A'^A. Assuming for instance that Chebyshev polynomials are used to approximate the 
inverse square root, the order of the Chebyshev polynomial is inversely pr^ortional to the square root of the smallest 
eigenvalue where the Chebyshev approximation has to work. In Table ^ we show the average smallest and eighth 
smallest eigenvalue of A^A with Dq being the Wilson Dirac operator, on different sets of gauge backgrounds with the 
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Figure 10: Flavor symmetry violation of staggered fermions on different gauge action backgrounds with different smearing 
transformations. All data on this plot is from quenched simulations with lattice spacing a ~ 0.17fm. 

same lattice spacing. The lattice size was 12'* and the (3 values were 6.0 (Wilson), 2.60 (Iwasaki), 8.40 (Symanzik) and 
1.04 (DBW2) and each set contained 50 independent configurations. For this qualitative test we fixed the value of s 
to be 0.5, except on the HYP smeared configurations, where it was set to zero. Although in principle s would have to 
be optimized for each type of gauge background separately, we chose to fix it close to the overall optimal value. This 
is sufficient for our purposes, moreover, an optimization for the smallest eigenvalues of A^A and for locality would 
yield different values. 

There is a clear trend that a more negative value of ci pushes up the smallest eigenvalues. However, the reader’s 
attention is called to the last entry in the table, where we have HYP-smeared the links in the fermion action and 
retained the Symanzik gauge action. The gain in time for simulating this action is almost a factor of two better than 
for the DBW2 action. 

Generally, overlap simulations are accelerated by projecting out (and treating exactly) the eigenvectors correspond¬ 
ing to the smallest few eigenvectors. This results in a gain of about a factor of 2 for the Wilson action, and smaller 
factors for the other actions studied. To facilitate a comparison we also included in the table the factor one can gain 
in speed compared to Wilson gauge action without projecting out any eigenvector. 

A major difference between the overlap and domain wall formulation is that for the overlap small eigenvalues of A'^A 
are only a nuisance, as they make the calculation more expensive. On the other hand, in domain wall simulations the 
extension of the lattice in the fifth direction is fixed and this results in different chirality violations and residual mass 
effects configuration by configuration. This is the main reason why improvement of the gauge action was so badly 
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Action 

(Ai) 

twijt 

(As) 

twi/t 

Wilson 

0.013(2) 

1.0 

0.061(1) 

2.17 

Symanzik 

0.044(4) 

1.84 

0.105(1) 

2.84 

Iwasaki 

0.065(4) 

2.24 

0.130(1) 

3.16 

DBW2 

0.160(5) 

3.51 

0.217(1) 

4.09 

Symanzik+HYP 

0.46(3) 

5.95 

0.737(2) 

7.52 


Table V: The average smallest eigenvalue (Ai) and the 8th smalles eigenvalue (As) of A and the factor of CPU time reduction 
compared to the Wilson action with no eigenvector projected out. The first two columns refer to the case when no eigenvalues 
are projected out, the last to columns to the one when the eight lowest eigenvectors are projected out and treated exactly. 
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Figure 11: The localization of the overlap in gauge backgrounds generated with various gauge actions. The parameter s was 
chosen to make the operator as local as possible in the given gauge background. 
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needed in the case of DW fermions. As our resnlts with the overlap snggest, improving the Domain Wall technology 
might be a better solution than going to extremes in tuning the gauge action. 


B. Localization 

We saw that the density of low modes of A^A depends on the gauge action. There have been specnlations on the 
connection between low modes of A'^ A and the locality of the overlap |^, It is therefore interesting to compare 
the locality of the overlap in the different gange backgrounds studied in the present work. This is done in Fig. |^, 
where the quantity 


\\Dovipix)\\ 


El 


f{r) = max 
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(27) 
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is plotted as a function of the so called “taxi driver distance”, r = \x^\ from a localized source ipkix) = 6(x)dkj. 

This quantity was introduced in |Q to measure the (non)-locality of the overlap operator. 

There are no surprises here. As seen in the figure, the general trend is that a large reduction in the number of 
small modes of A^A results in a slight improvement of the locality of the overlap, in accordance with the results of 
1^ , where different types and degrees of smearing were shown to have a similar effect. 


VIII. CONCLUSIONS 

In this paper we have investigated the reduction of lattice artifacts which can be achieved by altering the gauge 
field self-interaction and the fermion-gauge coupling. We focussed mostly - but not exclusively - on the reduction of 
lattice chiral symmetry breaking artifacts. Chiral improvement can be achieved by either altering the gauge action, or 
the fermion-gauge field coupling, or by a combination. That the two alterations produce similar effects is most starkly 
revealed by perturbation theory, where one sees that either effect alters the effective gluon propagator in fermionic 
Feynman diagrams. A variety of perturbative and nonperturbative tests reveal that a large ci in the gauge action can 
improve flavor symmetry violation for staggered fermions or the efficiency of implementing domain wall fermions or 
the overlap. However, a much more dramatic improvement can be achieved by replacing the thin link variable in the 
fermion action with smeared links. And the use of a large ci in the gauge action introduces a number of bad features 
into simulations: most notably a distortion of the heavy quark potential at short distance and long simulation time 
autocorrelations of the topological charge. 
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